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ABSTRACT 

An unexpected feature of Titan's atmosphere is the strong depletion in pri- 
mordial noble gases revealed by the Gas Chromatograph Mass Spectrometer 
aboard the Huygens probe during its descent on 2005 January 14. Although sev- 
eral plausible explanations have already been formulated, no definitive response 
to this issue has been yet found. Here, we investigate the possible sequestration 
of these noble gases in the liquid contained in lakes and wet terrains on Titan and 
the consequences for their atmospheric abundances. Considering the atmosphere 
and the liquid existing on the soil as a whole system, we compute the abundance 
of each noble gas relative to nitrogen. To do so, we make the assumption of ther- 
modynamic equilibrium between the liquid and the atmosphere, the abundances 
of the different constituents being determined via regular solution theory. We 
find that xenon's atmospheric depletion can be explained by its dissolution at 
ambient temperature in the liquid presumably present on Titan's soil. In the 
cases of argon and krypton, we find that the fractions incorporated in the liq- 
uid are negligible, implying that an alternative mechanism must be invoked to 
explain their atmospheric depletion. 
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Introduction 



A striking feature of the atmosphere of Titan is that no primordial noble gases other 
than argon were detected by the Gas Chromatograph Mass Spectrometer (GCMS) aboard 
the Huygens probe during its descent to Titan's surface in January 2005. The detected argon 
includes primordial 36 Ar present in subsolar abundance in Titan's atmosphere ( 36 Ar/ 14 N is 
found to be about six orders of magnitude lo wer than the solar va lue) and the radiogenic 
isotope 40 Ar, which is a decay product of 40 K ([Niemann et al.ll2005l ). The other primordial 
noble gases 38 Ar, Kr and Xe were not detected by the GCMS instrument, yielding upper 
limits of 10~ 8 for their atmospheric mole fractions. 

The interpretation of the noble gas deficiency measured in Titan's atmosphere has been 
the subject of several studies in the recent literature. Thus, lOsegovic &: Maxl (120051 ) pro- 
posed that these species could be preferentially stored in clathrates present on the satellite's 
surface. They calculated the c omposi t ion o f clathrates on the surface of Titan using the pro- 
gram CSMHYD developed by ISloanl ( 119981 ) and showed that suc h crystalline ice s t ructu res 
may act as a sink for Xe. However, the CSMHYD code used by lOsegovic fc Maxl (120051 ) is 
not suitable below 140 K for gas mixtures o f interest whereas the mean surface temperature 
of Titan is below 95 K (ICordier et al .1120091 ). and the authors did not explicitly calculate the 
trapping ef ficiencies of Ar and K r in cl athrates on the surface of the satellite. These consid- 
erations led iThomas et al.l ( 120071 . 120081 ) to rethink their results. In both studies, the authors 
found that the trapping efficiency of clathrates is high enough to significantly decrease the 
atmospheric concentrations of Xe and Kr irrespective of the initial gas phase composition, 
provided that these clathrates are abundant enough on the surface of Titan. In contrast, 
they found that Ar is poorly trapped in clathrates and that this mechanism alone could not 
explain the argon impoverishment measured in Titan's atmosphere. Another interpretation 
of the Ar, Kr and Xe deficiencies is that the haze present in Titan's atmosphere could simul- 
taneously tr ap these three noble gase s in a way consistent with the observed atmospheric 
abundances (IJacovi fc Bar- Nunl 120081 ) . In this mechanism, the open structure of the small 
aerosol particles would allow the noble gas atoms to fill their pores. All these hypotheses 
are based on different assumptions (requirement of large amounts of clathrates on the satel- 
lite's surface or formation of Titan's aerosols in exactly the same conditions as those used 
during laboratory experiments) that will need to be investigated by in situ measurements or 
observations performed by future spacecraft missions. 
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In this paper, we offer another hypothesis: that Titan's hydrocarbon lakes play a key role 
in the impoverishment of its atmospheric noble gases. Indeed, hundreds o f radar dark feature s 
interpret ed as hydrocarbon lak es have been detected in the polar regions ( IStofan et al.ll2007l ). 
Recently, ICordier et al.l (120091 ) have published a study of the chemical composition of Titan's 
lakes, which is based on the direct abundance measurements from the Gas Chromatograph 
Mass Spectrometer aboard the Huygens probe and recent photochemical models based on 
the vertical temperature pro file derived by the H uygens Atmospheric Structure Instrument. 
Here, we extend the model of lCordier et al.l (120091 ) by including simultaneously Ar, Kr and Xe 
in the composition of the liquid phase. We then explore the amount of liquid that is needed 
on the surface of Titan to account for the measured noble gas atmospheric abundances 
assuming solar abundances in the bulk system. 



2. Noble gases sequestration in Titan's surface liquid phase 



Here we consider the atmosphere and the liquid existing on the surface of Titan as a 
single, fully-coupled system. Our approach consists in i) computing the ratio iV NG /iV nitI . ogen 
where A^ng and A^ ni t r ogen are the total numbers of atoms of a given noble gas NG and of 
nitrogen, respectively, and ii) comparing this result to the ratio (N^Q/N nitrogen ) & derived from 
protosolar abundances (Lodders 2003). Our calculations always refer to nitrogen because it 
is the most abundant compound detected in Titan's atmosphere and because thermo dynamic 
equilibrium models predict that it is also present in the lakes (see lCordier et al.ll2009l ). Under 
these conditions, for a noble gas NG, we can write 



N- 



NG 



7V (liq) 4- N K 

iv tot,NG ' iv tot,NG 

Af(liq) _|_ AT-(atm) 

tot, nitrogen ' tot, nitrogen 



r(atm) 



iV, 



nitrogen 



where and iV^ 111 ^ are the total number of atoms of element NG or nitrogen in the 

liquid and in the atmosphere of Titan, respectively. 

The total mole number of element NG in Titan's liquid is given by 



iv tot,NG 



„(Hq) 

^NG 



M 



(2) 



where x^q is the mole fr a ction of atoms of NG in lakes computed following the method 



described in ICordier et al.l (|2009[ ). V^^ is the total volume of Titan's liquid in contact with 



the atmosphere, their mean density (in kg.m 3 ) and M the mean molecular weight of 
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the liquid. M is given by 



M 



Xj x Mj, 



(3) 



where the sum ^ . runs over all species present in the liquid phase. The composition of the 
liquid is calcu l ated v ia the thermodynamic equilibrium model of Titan's lakes described by 
Cordier et al.l ( 120091 ). In this approach, thermodynamic equilibrium, which translates into 
the equality of chemical potentials for each species listed in Table [TJ from N 2 to C2H 6 (except 
for H 2 ), can be expressed as 



vp,ii 



(4) 



where P is the total pressure at Titan's surface, yi and are respectively the mole fractions 
of the i compound in the atmosphere and the liquid, P vp , is its v apor pressure, and Tj i s 



its activity coefficient in the liquid determined with Equation 2 of iDubouloz et al.l ( 119891 ). 
Abundances of compounds below C 2 H 6 in Table [TJ are expressed proportionally to that of 
C2H 6 both in the precipitation and in the liquid existing in the lakes or in the putative 
porous network. 

In order to compute N^^ G for all noble gases, we have extended the thermodynamic 
equilibrium model b y adding Kr and Xe to the list of species already taken into account by 
Cordier et al.l ( 20091 ). The thermodynamic data essentially derive from the NIST databasqj 
in which the vapor pressures are expressed in the forms of Antoine's equations. On the 
other hand, the enthalpy of vaporization and molar volume of Ar, which are needed in 
the computation of its ac tivity coefficient via the determination of its solubil ity parameter 



(see iPrausnitz et al.lll986l ). have been updated relative to the value adopted bvlCordier et al. 



(2009) and now are derived from the laboratory measurements published by lFerreira fc Lobo 



( 120081 ) and iTegeler et al.l ( 119991 ). respectively. On the other hand, the enthalpies of vapor- 
ization and mola r volu mes of Kr and Xe all derive from the experimental data published by 
Ferreira fc Lobol (I2009I ). 



The total mole number of element NG in Titan's atmosphere is determined via the 
following vertical integration 



1 http : //webbook . nist . gov/chemistry 
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Table 1: Composition of liquid (mole fraction at given temperature). 





8 


7 K 




90 K 




93.65 K 






N 2 


1.22 


x 10" 


-2 


4.94 x 10" 


-3 


2.96 x 10" 


-3 




CH 4 


2.18 


x 10" 


-1 


9.74 x 10" 


-2 


5.56 x 10" 


-2 




Ar 


1.01 


x 10" 


-5 


4.95 x 10" 


-6 


3.09 x 10" 


-6 




Xe 


8.55 


x 10" 


-3 


1.52 x 10" 


-3 


3.09 x 10- 


-4 




Kr 


7.72 


x 10" 


-9 


3.13 x 10" 


-9 


1.92 x 10- 


-9 




CO 


1.24 


x 10" 


-6 


4.25 x 10" 


-7 


2.05 x 10- 


-7 




H 2 


2.94 


x 10- 


11 


4.08 x 10- 


11 


5.12 x 10- 


11 




C2H 6 


6.55 


x 10" 


-1 


7.62 x 10" 


-1 


7.95 x 10- 


-1 




CsH 8 


6.36 


x 10" 


-2 


7.40 x 10" 


-2 


7.71 x 10- 


-2 




C 4 H 8 


1.19 


x 10" 


-2 


1.39 x 10" 


-2 


1.45 x 10- 


-2 




HCN 


9.06 


x 10" 


-3 


2.08 x 10" 


-2 


2.89 x 10- 


-2 


(s) 


C4H10 


1.04 


x 10" 


-2 


1.21 x 10" 


-2 


1.26 x 10- 


-2 


(ns) 


C2H2 


9.83 


x 10" 


-3 


1.14 x 10" 


-2 


1.19 x 10- 


-2 


(ns) 


CH 3 CN 


8.48 


x 10" 


-4 


9.87 x 10" 


-4 


1.03 x 10- 


-3 


(ns) 


C0 2 


2.50 


x 10" 


-4 


2.92 x 10" 


-4 


3.04 x 10- 


-4 


(ns) 


CeH 6 


1.93 


x 10" 


-4 


2.24 x 10" 


-4 


2.34 x 10- 


-4 


(ns) 



Note. — From HCN to CeH6, compounds are in the solid state in precipitates and are assumed to dissolve 
when they reach the liquid phase, (s): saturated; (ns): non saturated. Ar is the total argon contained in all 
isotopes. 



(atm) 



tot.NG 



M( atm )(z) 



4tt (R 



Titan 



dz 



(5) 



where i?Titan is the radius of Titan, z the altitude and H the maximum elevation at which 
the Huygens GCMS started to collect data, the atmosphe ric density at the 



eleva tion z whose determination derives from the Huygens HASI data (jFulchignoni et al. 



2005|), ff^jr; ^ is the mole fra ction of element NG derived from the GCMS data at the ground 



level (INiemann et al.ll2005l ) and isassumed to be constant whatever the altitude. The mean 
molecular weight M^ tra \z) of the atmosphere is also derived from the Huygens GCMS data 
using an approach similar to the one used to calculate the mean molecular weight of the 
liquid (see. Eq. [3]). 
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Table 2: Assumed composition of Titan's atmosphere at the ground level. 



Atmosphere Mole fraction Determination 



H 2 


9.8 x IO" 4 


Huygens GCMS< a > 


N 2 


0.95 


This work 


CH 4 


0.0492 


Huygens GCMS^ 


CO 


4.70 x 10~ 5 


Cassim CIRS^ 


40 Ar 


4.32 x IO" 5 


Huygens GCMS^ 


36 Ar 


2.80 x IO" 7 


Huygens GCMS^ 


Kr 


io- 8 


Huygens GCMS® 


Xe 


io- 8 


Huygens GCMS^ 


C 2 H 6 


1.49 x IO" 5 


This work 



Note. — HOwen fc Niemann! |2009t) ; ( ^Niemann et al.1 (|2005h : Hde Kok et all (|2007h . N 2 and C 2 H 6 
abundances are determined from our model (see text). 



3. Assumptions about the Titan environment 

Atmospheric mole fraction data used as inputs in our model are gathered in Table [2j 
The atmospheric abundances of species which are not included in this table are assumed 
negligible. We have set the atmospheric mole fractions of Kr and Xe to 10~ 8 because they 
co rrespond to the dete ction limit of the Huygens GCMS instrument. Following the approach 
of ICordier et al.l ( 120091 ). atmospheric mole fractions of N 2 and C 2 H6 are treated as unknowns 
of our problem and atmospheric abundances of compounds below the one of C 2 H 6 are ne- 
glected. In order to investigate the temperature dependence of the mole fractions of noble 
gases trapped in Titan's surface liquid, we have allowed the ground temperature to range 
between 87 and 94 K. These extreme values bracket the mean ground tempera ture estimated 
to be ~90-91 K in northern polar regions where lake candidates are located (IJanssen et al. 



2009 ) and the tempe rature of 93.65 K measured by the Huygens probe at its landing site 
^Niemann et al-lbood ). 



Current inventories of Titan's lakes are estimated to range between 3 x 10 4 and 3 x 10 £ 



km 3 (ILorenz et al.l 120081 ) . In the present study, we consider a more generous range for the 



total volume of liquid existing at the surface of Titan, including a fraction that might reside 
in wet terrains. A substantial fraction of liquid incorporated in Titan's soil, even at the 
equator, is plausible because the Huygens probe firmly identified ethane in the mass spectra 
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taken from the surface (iNiemann et alJhood ). Moreover, it has been proposed that the 
subsurface of Titan is porous, implying that large amounts of li quid hydrocarbons coul d 
remain in contact with the atmosphere via open pores in the soil (IMousis fc Schmittl 120081 ) . 
Hence, assuming that the total volume of liquid existing on Titan's surface could be up to 10 
times the one estimated for the lakes, we consider here three different liquid volume values: 
3 x 10 4 , 3 x 10 5 and 3 x 10 6 km 3 . Note that the larger volume of liquid considered here remains 
~10 times lower than the total amount of ethane that could have precipitat ed on Titan's 



surfac e if one assumes a precipitation rate of 3.4 x 10 9 molecules. cm 2 .s 1 (ILavvas et al. 



2008al lbh over 4.5 gigayears 



10" 



^ 10 5 ' 



10" 



t — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 



Xe 




_i i i i i l i i i i i l i_ 



_i i l i i i i i l i i i i i l i i i i i l i i i i i_ 



87 88 89 90 91 

T(K) 



92 



93 



94 



4. Results 

Our calculations, including both 36 Ar and 40 Ar isotopes, produce a N\ T /-^nitrogen ratio 
(~2.24 x 10~ 5 ) remains more than three orders of magnitude lower than the value inferred 
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Fig. 1. — The global Nx e /N Q i tIogen ratio in Titan versus ground temperature for different 
volumes of lakes. The horizontal line represents the solar value. Other lines, from top to 
bottom, derive from our calculations for the considered volumes of liquid available at ground 
level. The assumed atmospheric mole fraction corresponds to the GCMS detection threshold 
(y Xc = 1(T 8 ). 

from solar abundances ((iVA r /iV nitrogen )Q = 5.26 x 10~ 2 ), irrespective of the adopted values 
for the t emperature an d total volume of lakes. Because 36 Ar is the main primordial argon 
isotope (Lodders 20031 ). we note that even if all the argon dissolved in the liquid was in the 



form of this isotope, this could not explain its apparent depletion in Titan's atmosphere. 
This behavior is explained by the very low solubility of argon in the C^Hg-dominated liquid 
in the considered temperature range. Even if the liquid volume is of the order of 3 x 10 6 
km 3 , almost all argon remains in the atmosphere, the fraction incorporated in liquid playing 
only a negligible role. Hence, whatever the source of argon, i.e. primordial or radiogenic, we 
conclude that the argon depletion observed in Titan's atmosphere cannot be explained by 
its sequestration in the liquid in contact with Titan's atmosphere. 

In the case of krypton, the computed A^r /-^nitrogen ratio (~5 x 10~ 9 ) is about 3-4 orders 
of magnitude lower than the solar value (( A^r /-^nitrogen )© = 2.83 x 10 -5 ), irrespective of 
the adopted values for the temperature and total volume of lakes. Similarly to the case of 
argon, dissolution of krypton in liquid hydrocarbon on the surface of Titan cannot explain 
the observed depletion. 

Figure [1] summarizes the behavior of the global ratio Nx c /N nitrogcn as a function of the 
temperature of Titan's surface and of the available volume of liquid. Because the thermo- 
dynamic equilibrium between the liquid phase and the atmosphere requires a larger mole 
fraction of xenon in the liquid, the resulting Nx e /N n i trogen ratio is larger than the one cal- 
culated from the detection limit of the GCMS instrument. Moreover, since the atmospheric 
mole fraction of xenon remains fixed in our system, its solubility diminishes in the liquid 
with an increase of temperature. 

Figured] shows that only the Nx e /-^nitrogen ratio reaches or exceeds the value determined 
from solar abundances in a temperature range consistent with the temperature inferred in 
polar regions and for a plausible amount of liquid on Titan's surface. Thus, a volume of liquid 
of about ~6 x 10 4 km 3 is required at 90 K to trap enough xenon to explain its apparent 
atmospheric depletion. The fact that xenon is much more soluble in liquid hydrocarbon 
than the other noble gases is consistent with their molar masses, i.e. ~39.95 g.mol -1 , 83.80 
g.mol -1 and 131.30 g.mol -1 for Ar, Kr and Xe, respectively. The most massive element 
remains the easiest to condense. 
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Discussion 



If one considers an atmospheric mole fraction of krypton lower than the detection limit 
(ukt = 10~ 8 ) of the GCMS instrument, the global ratio iVxr /^nitrogen should decrease linearly 
because N Kv / N nitvogcn — -^tot Kr /^nitrogen (the solubility of krypton appears to be very low: 



"tot,Kr 

f(atm) 



and N^ t ^ r oc i/kv (see Equation EJ) . In contrast, because xenon is much 



iV tot,Kr ^ iV tot,KrJ 

more soluble than the other noble gases in the liquid hydrocarbon, we have N^ Xe -A^xl- 
Assuming an atmospheric abundance of xenon lower than the GCMS detection threshold 
would also lead to a linear decrease of N% e /N nitrogeD because N^ Xe 3> N^l Xe anc ^ x x e oc yxe, 



implying that N^ Xe oc xx e (see Equation [2]). We find here a behavior close to Henry's law. 



The variation of pressure at ground level could also play a role in the dissolution of 
atmospheric compounds. However, Global Circulation Models (GCM) show that surface 
pressu re fluctuations due to wea ther conditions or Saturn's tidal effects should not exceed 
0.1% (ITokano fc Neubauerll2002T) . On the other hand, using the Cassini Synthetic Aperture 
Radar (SAR), IStiles et al.l (120091 ) have recorded surface height (referenced to the nominal 
2575 km radius) in the range -1500 m to +1000 m, yielding a maximum altitude difference 
of ~2500 m. We made a test at T = 90 K with a pressure 10% higher than the one measured 
by Huygens, i.e. a value corresponding to the unrealistic case where all the lakes and wet 
terrains are located at an altitude ~1500 m below the Huygens landing site. Even with this 
pressure variation, our calculations show that the fraction of dissolved noble gases is almost 
unchanged and that the afore-mentioned results remain unaltered. 

We conclude that noble gas trapping in Titan's hydrocarbon lakes and liquid contained 
in wet surfaces cannot be the unique answer to the problem of their apparent atmospheric 
depletion. The physical reality is probably composite and several combined effects might play 
a role with different efficiencies. For example, in order to explain the Ar deficiency, it has been 
proposed that Titan's building blocks were formed in a relative ly warm nebular environment 
which excluded both argon and molecular nit rogen (Owenlll982 ) or partly devolatilized du ring 
their migration within Saturn's subnebula (jAlibert &: Mousisl 120071 ; iMousis et al.l 120091 ) . If 
the lakes of Titan are the main sink of atmospheric xenon, then krypton must have remained 
sequestrated in the interior of Titan because none of the alternative trapping scenarios 
cited in the introduction predicts a krypton impoverishment relative to xenon. This study 
encourages direct measurements by future probes of the noble gas abundances in Saturn's 
atmosphere and/or in Titan's hydrocarbon lakes. 



We thank an anonymous reviewer for his constructive comments which helped us im- 
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